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We develop and analyze a theoretical model to study p-wave Feshbach resonances of identical 
fermions in atomic waveguides by extending the two-channel model of A.D. Lange et. al. [Phys. Rev. 

A 79, 013622 (2009)] and S. Saeidian et. al. [Phys. Rev. A 86, 062713 (2012)]. The experimentally 
known parameters of Feshbach resonances in free space are used as input of the model. We calculate 
the shifts and widths of p-wave magnetic Feshbach resonance of 40 K atoms emerging in harmonic 
waveguides as p-wave confinement induced resonance (CIR). Particularly, we show a possibility to 
control the width and shift of the p-wave confinement induced resonance by the trap frequency and 
the applied magnetic field which could be used in corresponding experiments. Our analysis also 
demonstrates the importance of the inclusion of the effective radius in the computational schemes 
for the description of the p-wave CIRs contrary to the case of s-wave CIRs where the influence of 
this term is negligible. 

PACS numbers: 32.60.+i,33.55.Be,32.10.Dk,33.80.Ps 


I. INTRODUCTION 

The recent progress in ultracold atomic physics 
provides exceptional possibilities for studying low¬ 
dimensional quantum systems. The confining geome¬ 
try of atomic traps can drastically alter their ultracold 
scattering behaviour near the so-called confinement in¬ 
duced resonances (CIRs) [lj - Employing magnetic Fesh¬ 
bach resonances [2], one can control the interaction be¬ 
tween the atoms and provide the conditions to experi¬ 
mentally observe the CIRs for identical bosons :3r[6j] and 
fermions (3-0, as well as distinguishable atoms [H|. CIRs 
have been extensively studied, e.g. in the context of 
bosonic s- flll - ITiii ] and d-wave [141] and dipolar [TS!, [3 
scattering, fermionic p-wave scattering mm , distin¬ 
guishable atom scattering [T3, fl8l - (2Cil ] in single-channel as 
well as in multichannel (transversal excitation) regimes 
[IH M H2 • Beyond this coupled 1-wave CIRs in cylindri- 
cally symmetric waveguides have been explored [ 23 , [24J . 

In spite of the impressive progress concerning the ex¬ 
perimental ||, fcj, !25| and theoretical [U, [261423 ] inves¬ 
tigations of CIRs the existing theoretical models need to 
be improved for a quantitative interpretation and guid¬ 
ing of the experiments. In the seminal work of Olshanii 
T] and in subsequent works jTl], 26, 27, 2f| 30], the sim- 


* saeidian@iasbs.ac.ir 

* melezhik@theor.jinr.ru 

4 pet er. schmelcher @physnet. uni- hamburg. de 


pie form of a pseudo-potential was used to model the in¬ 
teratomic interactions as compared to CIR investigations 
with more realistic interatomic potentials [ill , \l3i rul Ezt 
[ML (13, HU. However, all these approaches are of single¬ 
channel (internal) character, allowing one to explore only 
the main attribute of the Feshbach resonances in the 3D 
free space, i.e. the appearance of a singularity in the 
s-wave scattering length a s —> ±oo when a molecular 
bound state with energy Eb crosses the atom-atom scat¬ 
tering threshold at energy E = 0 in the entrance channel. 
However, other important parameters of the Feshbach 
resonance, such as the rotational and spin structure of 
the molecular bound state in the closed channel as well as 
the width A of the resonance characterizing the coupling 
T of the molecular state with the entrance channel, were 
ignored. In our recent work [32j the single-channel scalar 
interatomic interaction was replaced by the four-channel 
tensorial potential modeling resonances of broad, nar¬ 
row, and overlapping character for bosons, according to 
the two-channel parametrization suggested in [33]. This 
allowed us to calculate the shifts and widths of s-, d-, 
and g-wave magnetic Feshbach resonances of Cs atoms 
emerging in harmonic waveguides as CIRs. 

In the present work our particular focus are p-wave 
collisions of identical fermions due to their anisotropic 
character. Confinement induced shifts and widths of the 
resonances in an atomic waveguide could be envisaged, 
similar to what has been studied for s-wave interactions 
[32 ]. The main goal of the present work is therefore to 
extend our model [32] to the fermionic case. We then cal¬ 
culate the shifts and widths of the Feshbach resonance for 
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40 K atoms in its hyperfine state |F = 9/2, mjr = —7/2) 
for experimental conditions Q. The parameters which 
can be obtained from the experiments on magnetic Fes- 
lrbach resonances in free space [13, HH, namely, the spin 
characteristics, the background scattering volume Vb g , 
the resonant energies E c (or the corresponding resonant 
value of the strength B c of the external magnetic field), 
and the width of the resonance A (T) are used as input 
parameters of our model. 


II. FESHBACH RESONANCE MODEL IN FREE 
SPACE 

Let us first discuss an extension of our model [32|, de¬ 
veloped for resonant s-wave scattering, to the case of p- 
wave Feshbach resonances in 3D free space which can 
be observed in a single-component Fermi gas. Due to 
the Pauli exclusion principle the two-body wave function 
must be antisymmetric, and therefore only odd partial 
waves l contribute to the fermionic scattering process in 
case of a symmetric spin configuration. In contrast to the 
s-wave Feshbach resonance, the atoms near a p-wave Fes¬ 
hbach resonance have to tunnel through a centrifugal bar¬ 
rier to couple to the bound molecular state. Therefore, 
the wave function of the continuum can be significantly 
influenced by the bound state in the closed channel only 
in a narrow range of the magnetic field. Another feature 
of the p-wave Feshbach resonances is a doublet structure 
which arises from the magnetic dipole-dipole interaction 
between the electronic spins of the atoms. This splits the 
Feshbach resonance into distinct resonances indicated by 
their partial-wave projection m; = 0 and |to; | = 1 along 
the quantization axis. 

In this work we investigate collisions of two spin- 
polarized 40 K atoms in the hyperfine state \F = 
9/2, niF = —7/2). One can write the joint state of the 
atom pair as |/i = 9/2,771^ = —7/2)|/ 2 = 9/2, to/ 2 = 
—7/2 )|l = 1, mi) where mi = 0, ±1. The dipole-dipole in¬ 
teraction between the valence electrons of the atom pairs 
can be written as 

„ _ ^(r-SrXr •&)-(&■&) 

Uss — & ^3 

where a is the fine structure constant, Si is the spin of 
the valence electron of ith atom, r is the interatomic sep¬ 
aration, and r is the unit vector defining the direction of 
interatomic axis. It couples different partial wave compo¬ 
nents 1 1, mi) and \l',m[) with l' = l ± 2 or V = l ^ 0 [hit . 
and breaks the rotational invariance. The different com¬ 
ponents mi contribute differently in the dipole-dipole in¬ 
teraction, which means that molecular bound states with 
different mi possess different energies [34[. Consequently, 
the Feshbach resonances for different mi values couple to 
distinct molecular bound states and thus have different 
magnetic-field dependence. 

The difference in the pair collisions with the different 
mi projections can be understood intuitively by consider¬ 


ing the dipole-dipole interaction between the two atoms. 
The external magnetic field orients the spin of the va¬ 
lence electron of 40 K in |F = 9/2, tof = —7/2) state. 
When two dipoles are aligned side-by-side (head-to-tail), 
they are in a repulsive (attractive) configuration, corre¬ 
sponding to (r ■ Si) = 0 ((? • Si) = 1). 

The difference between the states to; = 0 and |to;| = 1 
is illustrated in Fig.3 of ref. [HI- In the case the motion 
is in a plane containing the magnetic field the interac¬ 
tion alternates between attractive and repulsive as the 
dipoles change between head to tail attraction and side 
by side repulsion. On the other hand, for the case that 
the atoms move in the plane perpendicular to the mag¬ 
netic held the dipoles are held in the side-by-side config¬ 
uration, and the interaction is only repulsive. Due to this 
effect one expects that the dipole-dipole interaction lifts 
the degeneracy between the to; = 0 and |m;| = 1 colli- 
sional channels which leads to a splitting of the Feshbach 
resonances, with corresponding resonance for |to;| = 1 at 
a higher energy (34}. This doublet structure disappears 
with increasing temperature due to a broadening of the 
resonance. 

For p-wave collisions the relevant parameter is the scat¬ 
tering volume, V p which has been parameterized as [T l 


k 3 cot5±(k, B) = 


1 k 2 
V P {B) + R(B) 


(1) 


near the magnetic Feshbach resonance, where V p and the 
effective range R were approximated as a quadratic func¬ 
tion of the magnetic field B 

-}— = c W+c^B + cMb 2 ( 2 ) 

V P\ B ) 


1 

W) 


d (0) + B + d {2) B 2 . 


( 3 ) 


The coefficients cW and dW were obtained by fitting to 
the experimental data and are given in Table I of |35| . 
A simple formula for an energy-dependent scattering vol¬ 
ume has been derived in the framework of multichannel 
quantum-defect theory (MQDT) 


V P (E, B) = Vb g (E) 


1 - 


A(l + #r) 




E 3 

eT 

*>g 


, ( 4 ) 


where Eb g 


h 2 

2 ^ a 't g 


is the energy associated with the 


background scattering length cib g = limF-u^IA^E)] 1 / 3 
(/i = to/ 2 is the reduced mass). Except for the case 
dbg A'* Rvdw , B b g ^ E V (i\\’ , and therefore E E vf iw 
(here R v dw = l/‘2(2fj,Ce/H 2 ) 1 ^ is the van der Waals ra¬ 
dius, C 6 is the corresponding van der Waals coefficient 
and E v dw = k 2 /2fj,R 2 dw ), when one does not need to 
take into account the effective range expansion of Vb g (E), 
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\mi\ 

B 0 [G] 

B*[G] 

$[MHz] 

A[G] 

t[^\\a bg [ao]\ 

0 

198.85 

178.508 

0.5684 

-20.342 

93.093 

-104.26 

1 

198.37 

175.414 

0.6950 

-22.956 

92.667 

-99.539 


the approximation (J7])-@. The width A of the magnetic 
Feshbach resonance reads 


.. r c~ 2 {E) . . 

lim —--— 7 — = —SuA 

e -¥0 2 tan Obg{E) 


( 10 ) 


TABLE I: The poles Bo and zeros B* of the scattering valume 
V P (B), the coupling strength T, the resonance width A, the 
magnetic moment difference and the background scatter¬ 
ing length dbg for a p-wave Feshbach resonance between two 
40 K atoms in the hyperfine state |F = 9/2, tuf = — 7/2) for 
two relative angular momentum projections mi along the axis 
of the magnetic field. 


Eq.fU) can be simplified to 


V p (E,B) = V b 


bg 


1 - 


B-B n - 


Sfi. _ 


(5) 


Fitting the expansion m to this formula one obtains the 
resonance parameters listed in Table I j35| . 

Using MQDT in the two channel case, the following 
formula for the phase shift in the open channel 


S(E, l ) = S bg (E , l) — arctan 




(E,l) 


E — E c + j tan A (E, l) 


_ ( 6 ) 
was obtained in ref.[35|. The first term 5b g (E,l) is the 
phase shift resulting from the scattering in the open chan¬ 
nel only, i.e. the background phase shift. The second 
term describes the resonant contribution originating from 
a bound state in the closed channel with the energy E c lo¬ 
cated close to the threshold of the open channel. C~ 2 (E) 
and tan A(-E’) are MQDT functions. The width of the res¬ 
onance r is multiplied by C~ 2 (E) which accounts for a 
proper threshold behaviour as k —> 0 . 

We assume that the bound state can be linearly tuned 
by a magnetic Zeeman shift, i.e. E C (B) = S/j,(B — B c ), 
where <5/r is the difference of the magnetic moments of 
the open and closed channels, and B c is the crossing field 
value of the bound state. In the case of a power law in¬ 
teraction potential r~ s the phase shift and MQDT func¬ 
tion exhibit the following Wigner threshold behavior as 

0 m 


and 


5 b g(E) 

T A t k 2l+1 

(7) 

2 and 



C~ 2 {E) - 

-T Bik 2l+1 

(8) 

tanA(-E) - 

tan A(0) 

(9) 


for all l. For van der Waals forces the higher order terms 
in k can be neglected in case E < E vd w- For alkali 
atoms E v dw ranges from O.lmK to 30mK [36}. Hence, 
in the ultracold regime (E < 1 fiK), one can safely use 


and the resonance position Bg, that is shifted from B c 
due to the coupling between the open and closed channels 
take on the following appearance 

r 

Bq = B c + —— lim tanA(-E). (11) 

26/1 E—>0 

Further, using the above parametrization, we extend 
the scheme suggested in ref. [33|] to describe the magnetic 
Feshbach resonances of the p-wave scattering in an ul¬ 
tracold 40 K gas. The two-body problem in free space 
permits the separation of the center-of-mass and relative 
motion yielding the following Hamiltonian for the relative 
atomic motion 


H(r,9) 



I + V{r) 


( 12 ) 


Here V (r) is the two-channel interatomic potential and 
I is the unit matrix and 7' is the relative radial coor¬ 
dinate. Let us suppose that initially the spin-polarized 
atoms are prepared in the entrance channel |e) and the 
“closed channels” |c) supports molecular bound states at 
Bo. The quantum state of an atomic pair with energy E 
is then described as \ip) = ip c (r)\c } + tp e ( r)|e) satisfying 
the Schrodinger equation with the Hamiltonian (fl2l) . A 
two-channel square-well potential 


“) « r <“> < 13 > 


( 00 0 \ ,. , 

= { 0 0 jW r>fl ) 

is employed to describe the colliding atoms in the “en¬ 
trance channel” |e) and the weakly-bound molecule in the 
“closed channel” |c) near a Feshbach resonance. The cou¬ 
pling constant hfl induces Feshbach couplings between 
the channels. For r < a, we assume that the attrac¬ 
tive potential supports multiple molecular states - that 
is 14,14 > E vdw = h 2 j2[iR 2 vdw ( R vd w = 65 a 0 and a 0 
is the Bohr radius). For r > a, entrance- and closed- 
channel thresholds are set to be E = 0 and E = 00, re¬ 
spectively. Here a = (U) 1/3 and V = R ldw is 

the mean scattering volume f37j , and T(a:) is the gamma 
function. For the pair of 40 K atoms we have the value 
a = 63.4ao- 

Such a choice of the interatomic interaction permits 
a simple parametrization of the atom-atom scattering in 
universal terms of the energy of the bare bound state E c , 
the Feshbach coupling strength T of the bound molecu¬ 
lar state with the entrance channel and the background 
scattering volume V b g , which is convenient for an analysis 
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of experimental data near magnetic Feslibach resonances 

a fl. 

When the mixing between the closed channel and the 
entrance channel is weak and the background scattering 
length |a{, g | considerably exceeds the range of the inter¬ 
atomic interaction a, the p-wave scattering volume V p 
obeys the following expression 


1 _ 1 F /2 

V p -V ~ V bg - V + VE C 


(14) 


The parameters Vb g , 5/j,, Bq, and A from jUJ together 
with V P {B) given by Eq.©) are used for fitting the diag¬ 
onal terms V c and V e in the tensor potential (1131) . The 
nondiagonal terms Htt are defined by T/2 = 2 6 2 V Cl where 
tan 29 = 2hfl/(V e - V c ). 

The scattering volume V P {B) is then calculated for dif¬ 
ferent B and varying parameters of the potential V by 
solving the Schrodinger equation 


2 /j 


I + B + V(r) } \i/;) = E\i/>) (15) 


with scattering boundary condition 

Mr) -»• e ifez + e ikr ,Mr) ->• 0 (16) 

r 

at kr —» oo for the fixed E — > 0 (k = \J2[lE/h) [Hj]. 
The diagonal matrix B in Eg. (1 15 [) is defined as B cc = 
5[i(B — B c ) and B ee = 0. After separation of the angular 
part in Eq. m we come to the system of two coupled 
radial equations 


2/i dr 2 


H 2 l(l + 1 ) 

2 fir 2 


Bn 


M r )+Yl 

(17) 


with l = 1, and i,j = e, c. 

By Settin § tan 45) = ~ and Varying Vc ’ V * 

and Q we obtain an excellent agreement of the calculated 
p-wave energy dependent scattering volume V P (B,E) = 
— k~ 3 tanSi(B,k) (Eq.©) with the expression © from 
[ 35 } for 40 K atom in the hyperfine state | F = 9/2 : itif = 
-7/2) for 196G < B < 202G. The found zeros B* 
and poles Bq of the scattering volume and the coupling 
strength T are given in Table I. Here B c = Bq + /3A, 
A = B* — B 0: and /3 = Vb g /Vb g — V. We note that 
such a procedure yields a fc 3 cot Si(B , k) coinciding with 
the experimental data [34] for low colliding energies at 
T ~ 1—10 nK. The agreement, however, weakens with in¬ 
creasing k. Therefore, our 1/R = -l/2-^ 7 (l/V p (B,E)) 
differs from the value given in (34j by about a factor of 
3. This fact is, however, understandable because we did 
not use 1/R directly from [ 34 } hr our fitting procedure 
for obtaining the parameters of the interaction potential 

m 

Fig© shows our results obtained for the state mi = 0 
for the scattering volume as a function of the magnetic 


field B at the temperature T = 1.0 nK. For compari¬ 
son we have also plotted the analytical curve (5) which 
is in a good agreement with the numerical result. The 
computations were performed for V e = 1.5 x 10~ 2 [a.u.], 
V c = 5.9 x 10 _ 2 [a.u.] and fl = 1.2 x 10 _ 6 [a.zi.] modeling 
unbound two-atonr states in the entrance channel |e > 
and the resonance state with the parameters fixed for 
mi = 0 in the Table I. 


III. FESHBACH RESONANCE MODEL IN A 
WAVEGUIDE 


Next we analyze the scattering properties of the p- 
wave magnetic Feshbach resonances in harmonic waveg¬ 
uides permitting unbound motion in the longitudinal z- 
direction and a transversally strongly confined p-motion 
in the potential 1/2pw[j_p 2 . To describe the scattering 
process in the waveguide we have to calculate the scat¬ 
tering amplitude f P {E) by integrating the Schrodinger 
equation 


H 2 _ 2 1 

2 t V ‘ + 


2 „2 


/ + B + V(r) 


with the scattering boundary conditions 


m 

(18) 


M r ) = (sm(k 0 z) + sgn(z)f p e lk °^ $ 0 (p) , Pc,i(r) 0 

(19) 

at | z |=| rcosB |—> 00 adopted for a confining trap 
0 - Here, $ 0 (p) is the ground-state wave-function 
of the two-dimensional harmonic oscillator and ko = 
\j2[i{E — huj±)/h = ^/2iiE\\/h. We consider pair col¬ 
lisions of identical fermionic potassium atoms, therefore 
j(jfc)ie scattering wave-function | ip) must be antisymmetric 
with respect to the exchange 2 —> — z . In the presence of 
the harmonic trap 1/2pcu^p 2 the azimuthal angular part 
of the solution is separated, and Eq. (fT5|l is reduced to the 
coupled system of two 2D Schrodinger-type equations in 
the plane {p, z}. To integrate this coupled channel scat¬ 
tering problem in the plane {r, 0} we have extended the 
computational scheme developed in [13]. The integra¬ 
tion was performed in the units of the problem leading 
to the scale transformation: r —> =, E ~>■ J/, F — > 

and uj± -A ^ with E 0 = (= 6.8 x lO" 9 ^ for 40 K 

atoms), and ojq = ^ (= 2.8 x 10 5 kHz for 40 K atoms). 


IV. COMPUTATIONAL APPROACH 

First, we discretize the Schrodinger Eq. (fTS]l on a grid of 
the angular variable 1 , which has been introduced 

in 0. We expand the solution of Eq. GHD in the basis 
9j{0) = X0V 1 Pi{cosd)(P~ 1 )i J according to 

Nq 

\P) = V’iM) = ~Y1 Y^9j{0)u\ l) (r) , (20) 

i=e,c i=e,c j =1 
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where P -1 is the inverse of the Ng x Ng matrix P with 
elements defined as Pji = XjPi(cos6j) and A j are the 
weights of the Gauss quadrature. Substituting dm) into 
(H51l results in a system of 2 Ng Schrodinger-like coupled 
equations with respect to the 2 x Ng dimensional un¬ 
known vector u(r) = {uji(r)} = {X) /2 u!f {r)}^ 1}i=e>c 

\K[x) + 2 (El - U{r{x)) - B)]u(r(x)) = 0 (21) 

with 

Kf;.[x) = -7^)- 

N e -l 

-Y'PjiKi + iKP-'MSu' ( 22 ) 

z=o 

and 

u jj'( r ) = [Vn'( r ) + PjSw]Sjj' ,Pj=r sin Oj . (23) 

with j3{x) = (e 7 — l)/(r m 7 e 7X ) and we mapped and dis¬ 
cretized the initial variable r £ ( 0 , r m ] onto the uniform 
grid Xj £ ( 0 , 1 ] according to 

e 7*j _ i 

r j = r m eT _ ) 3 ~ h 2) •••) N , (24) 

where r m is chosen in the asymptotic region r m —> oo 
and 7 > 0 is a tuning parameter. Using the finite dif¬ 
ference approximation, the above Schrodinger-like equa¬ 
tions are reduced to a system of algebraic equations ac¬ 
cording to ref. [u|. By solving this system of equations 
for fixed colliding energy E and matching the calculated 
vector il> e (E,r = r m ,9j) with the asymptotic behavior 
m at r = r m , we find the scattering amplitude f p . 
from which one can obtain the observable transmission 
coefficient T{E) = |1 + f p (E )| 2 in the waveguide. 

V. RESULTS AND DISCUSSION 
A. Transmission coefficient 

We have performed two sets of integration of Eq. da 
for varying interatomic interaction (defined by a vary¬ 
ing magnetic field B) of 40 /\ atoms in the hyperfine 
state |F = 9/2, mp = —7/2) and the relative angu¬ 
lar momentum state 1 1 = 1, to; = 0) near the Fes- 
hbach resonance Bq = 198.85G: the dependence of 
the transmission T(B,lo j_,e) = |1 + f p (B,w±, e )| 2 on 
the trap frequency w_l, defining the transversal volume 
V± = _l)) 3 / 2 , and the rescaled longitudinal energy 

e = 2 K£l = = itr ~ k was analyzed. Calcu¬ 

lations were performed with the parameters V e , V c , and 
Q of the interaction potential V (11311 fixed in Section II 
by fitting the parameters of the Feshbach resonance of 
potassium atoms in free space. 



FIG. 1: (Color online) The p-wave scattering volume V v cal¬ 
culated for potassium atoms near the Feshbach resonance in 
the hyperfine state |F = 9/2, m_F = —7/2) and the relative 
angular momentum state |Z = 1, mi = 0) as a function of the 
magnetic field B. The solid curve shows the analytical result 
(5) and the dots show the numerical result (see text). 




FIG. 2: (Color online) The transmission T in the harmonic 
waveguide with (a) o>_l = 0.002 and (b) oj± = 0.006 as a 
function of external magnetic field B for several rescaled lon¬ 
gitudinal energies e. 


In Fig (2] the dependence of the transmission 
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FIG. 3: (Color online) Transmission T in the waveguide for 
the low energy limit e = 5x 10 -4 as a function of V p /Vo_. 



FIG. 4: (Color online) The position of the p-wave CIR in the 
harmonic waveguide as a function of e predicted by Eqs. (ESI) 
(solid curve) and (1251) (dashed curve). 


T(B,ui±,e) on the external magnetic field B is given 
for two different u>± and varying colliding energy. The 
p-wave Feshbach resonance is manifested in the waveg¬ 
uide as a minimum in the transmission T(Bcir) = 0 
(i.e. as a p-wave cir) E! mm. Fign demonstrates 
the strong dependence of the shift of the resonance 
Bo — Bqir in the confining trap on the trap frequency 
wj_ similar to the s-wave CIRs |32|]. This effect, as in 
the case of the s-wave CIR II , can be used to control 
the p-wave CIR position. However, for a complete 
control of the p-wave resonance shift in the waveguide 
one has also to understand the strong dependence on 
the energy e shown in Fig(2] This effect has to be 



FIG. 5: (Color online) The position of the p-wave CIR as a 
function of e in harmonic waveguides with cj_l = 0.002 and 
0.02. The curves show the analytical results (|29) and the dots 
correspond to the numerical results. 



FIG. 6: (Color online) The position of the p-wave CIR in the 
harmonic waveguide as a function of u}±. The curves show 
the analytical results EH) and the dots correspond to the 
numerical results. 


contrasted to the position of the s-wave CIR given by 
a s (BciR)/a± = 0.68 (where a± = and a s is the 
s-wave scattering length) [lj which is taken in the e —> 0 
limit and was confirmed in recent experiments H. The 
strong dependence of the p-wave CIR position on e is 
equally known 0- However, subsequent investigations, 
both of analytical ill, 29] and numerical character 13 . 
have demonstrated strong deviations from the results 
in ref. 0- In the next subsection we shall discuss 
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FIG. 7: (Color online) The width of the p-wave CIR of potas¬ 
sium atoms in the hyperfine state |F = 9/2,m_F = —7/2} 
and the relative angular momentum state \l = 1 , mi = 0) as 
a function of the rescaled energy e in harmonic waveguides 
with wj_ = 0.002 and 0.02. N = 100 and 1 for u>± = 0.002 
and 0.02 respectively. The curves show the analytical results 
m and the dots show the numerical results. 



FIG. 8: ((Color online) The width of the p-wave CIR of potas¬ 
sium atoms confined by harmonic waveguide in the hyperfine 
state |F = 9/2, tuf = —7/2) and the relative angular mo¬ 
mentum state \l — 1, mi = 0) as a function of the waveguide 
frequency wi for e = 5.0 x 10~ 4 and 9.0 x 10 _1 . N = 100 
and 1 for e = 5.0 x 10” 4 and 9.0 x 10 _1 respectively. The 
curves show the analytical results m and the dots show the 
numerical results. 



FIG. 9: (Color online) The real part of the scattering am¬ 
plitude f p (dots) as a function of the V p /V± along with the 
analytical results (1301 (curves) for e = 0.9. 

the origin of this difference. To this end we study the 
transmission T(V P (B), wj_, e) =| 1 + f p {V p (B), u±, e) | 2 
as a function of the p-wave scattering volume V P (B)/V± 
“normalized” to the transversal volume V± and as a 
function of the trap fequency in the low energy limit, 
which is shown in Fig|3] The dependence of the p-wave 
CIR position Vf IR /V ± , defined as T(V p CIR /V ± ) = 0, 
on ui _l qualitatively shows the same behaviour as in our 
previous work 13] . However, the calculated values of 
K cm /v ± | in [131 ] are approximately two times larger 
than the present values. This difference stems from the 
definition of V p for the chosen interatomic interaction 
in [l3}. Specifically, in ]13] the effective range R was 
not included in the definition of the scattering volume 
V p and the interatomic interaction was of single channel 
character, contrary to the tensorial structure of V in 
m- We note that we have shown recently (HI using a 
tensorial potential of type m that both these effects 
do not influence the position of the s-wave CIR. 

B. Shift of p-wave Feshbach resonance in harmonic 
waveguide 

In 0 the CIR of a pair of identical fermions has been 
studied analytically by means of a fermion-boson map¬ 
ping. In this way the authors have obtained the ratio 

V CIR 1 

^ = [12C(--,l-e)]- 1 (25) 

for the position V p IR of the CIR corresponding to the 
minimum in the transmission T. This ratio depends only 
on the rescaled longitudinal energy e and the dependence 
is defined by the Hurwitz zeta-function £. In Fig0] we 




























show the resonance position predicted by this formula as 
a function of e (the solid curve). On the other hand, in 
a the dependence of the p-wave CIR position on e is 
given by the simple expression 

yCIR 

-^ 7 — = -[4(1 - e )]- 3/2 , (26) 


predicting a smoother dependence on e and a strong de¬ 
viation from (12511 with varying e (see the dashed curve in 
Fig©). Thus, both formulae exhibit a strong dependence 
of the p-wave CIR position on e in support of our above 
result (see Fig© but do not describe the strong depen¬ 
dence of the ratio V p IR /V±_ on the trap frequency uj± 
we find: see Fig© and Figs.lO(a,b) in our previous work 
[Dl . Our above analysis has shown that the origin of the 

drawback of both formulae m and JUD is the neglect- 

■ 2 

ing of the second term in Eq. (JT|) in the derivation 
of these formulae. 

In a recent work [39 ] it was concluded that, unlike in 
the case of s-wave interaction, the p-wave effective range 
is essential in the strongly interacting regime i.e. near the 
corresponding Feshbach resonance and the shift of the 
p-wave magnetic Feshbach resonance in confining har¬ 
monic waveguide was calculated analytically. Assuming 
that the external magnetic field lies in the ccz-plane with 
angle 4> with respect to the 2 -axis, and that the incom¬ 
ing particles are in the ground state of an axially sym¬ 
metric harmonic confinement, the authors of 1391 ] have 
obtained an analytical expression for the scattering wave 
function. They arrive at the energy dependent ID scat¬ 
tering length: 

tan<5i 6a±(D{® cos 2 </> + sin 2 (f )) 

ko Dl°D\v cos 2 (j) + D^Dl? sin 2 

(27) 

where k 0 = sj2/i{E - fiw±)/h = ^ (e = - ±), 

= D™ 12C(-|,1 - e), D™ = D™ + C* D (e), 

D^ and C 4D (e) is a function defined in 

A 

For anisotropic interaction (Dq D ^ D\ D ) with <j> = 0 
or (f> = j as well as for isotropic interaction (Dg D = 
D{ D ), one obtains: 


O-ID 


6u_l 

D 1D 

mz 


(28) 


which yields a resonance (Z)^ = 0 ) at — 12 £(— 1 — 

e) = 0 , from which one obtains for the resonance condi¬ 
tion i.e. for the position of the p-wave CIR: 



[^ + 12C(-i 1-e)]- 1 


(29) 


This formula is essentially different from Eas. E5l) and 
(EH 1 ) by the presence of the term . 

In our case of two spin-polarized 40 K atoms | F = 
9/2, niF = — 7/2) with relative quantum numbers l = 1 


and m; = 0 trapped in the waveguide and equipped with 
a longitudinal magnetic field we expect the resonance po¬ 
sition to obey Ea. (l29l) . The results presented in Fig© 
demonstrate excellent agreement of our calculation of the 
p-wave CIR position V p IR /V±_ with the analytic result 
(EH 1 ) in the complete region below the first threshold of 
transverse excitation e = 1 , which persists in a broad 
range of varying io±. Note also the emerging stronger 
than linear dependence on the energy e with increasing 

Wj_. 

Fig© shows the resonance position as a function of the 
trap frequency to± for the zero energy limit e = 5 x ICC 4 
as well as near the first threshold of transverse excitation 
e = 0.9. It demonstrates the close to linear dependence of 
the resonance position V^ IR /V± on v / u7T as w_l —> 0 re¬ 
gardless of the value of the energy (see figure inside). The 
numerical results are in good agreement with Eg. (1291) . 


C. Width of p-wave Feshbach resonance 


Next we analyze the widths of the p-wave CIRs in 
harmonic waveguides. Since at the resonant field Bqir 
(corresponding to the position of CIR) the transmission 
reaches T(Bcir) = 0 and the maximal T value is 1 we 
define the width Aip = B + — R_, where B + and are 
the fields correspondingly right and left to Bcir and the 
transmission approaches the value T(B±) = 1/2. This 
definition differs from that for a free space magnetic Fes¬ 
hbach resonance where it reads B* — B 0 [33i]. By using 
the analytical formula for the scattering amplitude in ref. 


f P {B) = 


_- ik 0 _ 

a i d(B) + i?fcg + iKo 


(30) 


valid near the CIR in the zero-energy limit kg = 
y/2/ZE|i"/7i = y/2fi(E — —> 0 and Ea. (l28l) we ob¬ 

tain 


Ai d = A] 


a+ 

ct+ — 1 



(31) 


where a± = 6 Q ^ a ^ g ( g ±) + lycra- In Figs. □ and © 
we show the calculated resonance width Aid as a func¬ 
tion of the rescaled energy e and waveguide frequency 
w_l, respectively, together with the analytical result ac¬ 
cording to Eq. m ■ One observes a good agreement for 
uj± < 0.01. However, we encounter major deviations with 
increasing co± except for the zero energy limit e —> 0 . 
These deviations are due to the fact that the analytical 
formula (1501) which we used to derive (1311) , was obtained 
in the zero energy limit and does not work for larger en¬ 
ergies. In Fig© we present the real part of our scattering 
amplitude Re{f p ) along with the analytical results from 
as a function of V//Vj_ for e = 0.9 (large energy) 
for comparison. As the trap frequency w_l increases we 
observe deviations between the analytical and numerical 
results for f p . 


















9 


The same way as in the case of the free space reso¬ 
nance, the width Ai d of the CIR narrows with decreas¬ 
ing energy e (see Fig[7]). Fig[8] shows that in a harmonic 
waveguide there is a region where we have a possibility 
for narrowing the width by decreasing the trap frequency. 

VI. CONCLUSION 

We develop and analyze a theoretical model to study 
Feshbach resonances of identical fermions in atomic 
waveguides by extending the two-channel model sug¬ 
gested in |33] and adopted in Q for confined bosons. 
In this model, the experimentally known parameters of 
Feshbach resonances in free space are used as an input. 
Within this approach we have calculated the shifts and 
widths of p-wave magnetic Feshbach resonance of 40 K 
atoms in the hyperfine state |F = 9/2, mp = —7/2) and 
for the relative angular momentum state |Z = l,mj = 0) 
emerging in harmonic waveguides as p-wave CIRs. We 
find a linear dependence of the resonance position on the 
longitudinal colliding energy below the threshold for the 
first transverse excitation. It is shown that in a har¬ 
monic waveguide there is the possibility to decrease the 
width of the p-wave Feshbach resonance by decreasing 


the (transversal) trap frequency which could be used in 
corresponding experiments. Our analysis demonstrates 
the importance of including the effective range terms in 
the computational schemes for the description of the p- 
wave CIRs contrary to the case of s-wave CIRs where the 
impact of the effective radius is negligible. In previous in¬ 
vestigations of the p-wave CIRs in harmonic waveguides 
0, S3, 0 the effects due to the effective range have 
been neglected. The developed model can be applied for 
a quantitative analysis of other p-wave CIRs following a 
different spin structure and for confining traps of differ¬ 
ent geometry including effects due to anharmonicity and 
anisotropy. 
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